February 1, 2008 2:18 WSPC/INSTRUCTION FILE qmposter 



International Journal of Modern Physics E 
© World Scientific Publishing Company 



PRODUCTION OF ANTICENTAURO EVENTS IN 
ULTRA-RELATIVISTIC HEAVY ION COLLISIONS 



G. Sood 

Physics Department, Government College (GCG -11), 
Chandigarh, 160011, India, 
gopika&hepphys . com 

Received (received date) 
Revised (revised date) 

We propose a novel method for studying the production of anticentauro events in high 
energy heavy ion collisions utilizing Chebyshev expansion coefficients. These coefficients 
have proved to be very efficient in investigating the pattern of fluctuations in neutral 
pion fraction. For the anticentauro like events, the magnitude of first few coefficients is 
strongly enhanced (~3 times) as compared to those of normal HIJING events. Various 
characteristics of Chebyshev coefficients are studied in detail and the probability of 
formation of exotic events is calculated from the simulated events. 



1. Introduction 

In ultra-relativistic heavy ion collisions the rapid expansion of collision debris leads 
to the production of vacuum states with anomalous chiral order parameters. As a 
result of this, interior vacuum straighten out and radiate away its pionic orientation. 
If the deflection of vacuum orientation is in the tt'^ direction all the condensate 
radiation will be tt". On the other side, if deflection is orthogonal to the 7r° direction 
all the emitted pions will be charged. The former case is called "Anti-centauro" 
behavior and the latter is "Centauro" . Both of these type of events are detected in 
the cosmic ray experiments^. Various models^ have been developed to understand 
the characteristics of these events and numerous signatures ^ for identification of 
such events are suggested. 

The LHC facility SI at CERN will soon be operational and it would be interesting 
to search for anomalous events of centauro/anticentauro nature at LHC energies. 
Further, it is expected that at these high energies the production of exotic events 
will be enormous. So, it is useful to perform a theoretical simulation to identify 
the specific variables and mathematical techniques that may be informative in this 
study. The primary and widely accepted signal for observing these events is to 
look for anomalous fluctuations in neutral pion fraction (correspondingly, large 
fluctuations in the energy ratio, electromagnetic to hadronic). In this paper, we 
introduce a technique of Chebyshev coefficients ^ to study the fluctuations in the 
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neutral pion fraction. 

2. Characteristics of Chebyshev polynomials 

Chebyshev polynomials are actively used in the field of mathematics for computing 
the derivative/integral of a non-singular function in the finite domain and in solving 
various differential equations In the High Energy Physics area, this method is 
applied in the track reconstruction procedure^, vertex calculation and in studying 
the pattern of fluctuations The advantage of Chebyshev expansion technique 
is that it can be extended to any variable and it can be easily applied to any 
experimental data. In this paper, we have used this novel technique to study the 
fluctuations in the azimuthal distribution of neutral pions to the total number of 
pions produced in an event. The neutral pion fraction, /(</>), can be expanded in 
terms of sum of Chebyshev polynomials as : 

C ^ 

f{'l^) = ^ + Y.^M<P) (1) 

fc=l 

where Tfc((/>) is the Chebyshev polynomial of degree k. The coefficients Cj, j=0, 1, 
2, ■ ■ ■ , N-1, of Chebyshev expansion are defined as 

2 ^ 
fc=i 

^i = N^ N ^^''"''^ N ^ ^ ^ 

fe=i 

The coefficients of Chebyshev polynomial carry all the information on the fluctua- 
tion of The useful measures are provided by ranking the coefhcients, Ck, in 
order of the magnitude of the absolute values, i.e., 

I Co I > I Ci I > I ^2 I > I C3 I • • • (4) 

The magnitude and range of these coefficients define the strength and level of 
fluctuations present in the data sample. 



3. Simulation 

The cosmic ray evidence on "anti-centauro events" serve as a great motivation for 
performing an exotic event search at Large Hadron Collider (LHC). Therefore, for 
any study, in this direction, it is essential to follow the hints and inferences provided 
by the cosmic ray experiments. We have taken the spectacular JACEE prototype 
^ as a reference for generating the anticentauro events artificially. The standard 
Monte Carlo program, HIJING, is used to simulate minimum biased (0-20 fm) Pb- 
Pb collisions at LHC energies. Default values are taken for all the parameters. Two 
ensembles of 10 K events, (i) the normal HIJING events and (ii) the HIJING events 
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Fig. 1. First Chebyshev coefficient, Co, for HIJING (dashed line) and anticentauro like events 
(solid line). 



with anticentauro type fluctuations embedded in it, are generated. The anticentauro 
event is modeled according to the {l/2y/J) distribution. For the present simulation, 
a domain of size 60° is injected into the normal HIJING event. This is established 
by flipping the charged pions into neutral pions unless they attain a value of the 
neutral pion fraction, / > 0.9. 



4. Analysis 

We calculate the neutral pion fraction, /, for the generated data sets (HIJING as 
well as "anticentauro like events" ) in the Arj — A<p phase space as, 

/= (5) 

This calculation is performed on an event-by-event basis for various azimuthal bins, 
viz, A(j) = 27r/18, 2tt/9, 27r/6, 27r/4, 2tt/3, tt and 2tt and in the broad pseudorapidity 
region (77 = -10 to 10). The Chebyshev coefficients corresponding to the computed 
neutral pion fraction are determined according to the equation 3, as described in 
section II. Various characteristics of Chebyshev polynomials are investigated in view 
of their ability to pick up the events with anomalous fluctuations in the neutral pion 
fraction. The details are described below : 



4.1. First order coefficients 

The first Chebyshev coefficient, Co, for = 27r/18, for normal events (/ = 0.33) is 
plotted in Fig. 1. This distribution is Gaussian (dotted line) and it is chosen as ref- 
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Table 1. Mean value of the first coefficient, Cq, for the exotic 
and normal HIJING events. 



Mean Co 

Azimuthal size Anticentauro like event HIJING event 



27r/18 


124.8 


45.0 


27r/9 


63.2 


23.0 


27r/6 


42.7 


15.4 


27r/4 


28.5 


10.2 


27r/3 


21.4 


7.7 


TT 


14.3 


5.2 


27r 


7.1 


2.6 



erence for studying the pattern of fluctuations in the generated anti-centauro type 
events. This is an ideal choice for comparison as there are no hidden fluctuations 
present in this sample and it purely represents the true nature of the system. In the 

same figure, wc have plotted the first order coefficient (Co), for the anticentauro 
like events (solid line). The absolute value of the magnitude of this coefficient is 
large as compared to that of normal events. Similar behavior is observed for other 
order coefficients, i.e., for the exotic events the value of coefficients is significantly 
enhanced as compared to the normal events. Table 1. describes the detailed results 
of the variation of absolute magnitude of first order coefficient, Co, with azimuth, 
for the top 5% central events. The mean value of the coefficient, Cq, for anticen- 
turo like events is approximately thrice the mean value of normal HIJING events. 
It is observed that with increase in size of domain the magnitude of coefficients 
decreases. This behavior is seen for both HIJING and anticentauro like events but 
the magnitude of the coefficients for exotic events is much more than normal events 
for all the azimuthal bins. 



4.2. Higher order coefficients 

The coefficients of higher degree Chebyhshev polynomial, for 0< k < N-1 contain 
all the information of fluctuations. In our study, we have calculated ICQ coefficients, 
Co to C99. All these coefficients, Cfc (k = 0, 99) are pure Gaussian (for normal 
HIJING events), however, the absolute values of their magnitude decreases very 
rapidly. Fig. 2 shows the distribution of some higher order coefficients, C30, C40, C50 
and C90. It is observed that higher order coefficients gradually attain the null value. 
Therefore, the Chebyshev expansion scries can be truncated at some order m < N 
and the extreme higher order coefficients are neglected. It is important to note that 
the absolute value of the magnitude of higher order coefficients for the "anticentauro 
like events" is still much higher than the corresponding HIJING events as observed 
for the flrst order coefficient. 
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4.3. Correlation of Chebyshev coefficients 

The magnitude and range of the Chebyshev coefficients is very useful in predicting 
whether or not the fluctuations in the neutral pion fraction are of non-statistical 
nature. A more critical criteria can be evolved by considering the two coefficients 
simultaneously. For this, it is necessary to study the correlation between various 
coefficients. Fig 3. shows the relationship between the two consecutive coefficients 
Ci and C2. It is seen that the coefficients are correlated i.e., with increase in Ci, 
the coefficient C2 also increases. Similar behavior is observed for other coefScients. 
The corresponding magnitude of coefScients for HIJING events is very less. 
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Fig. 2. The distribution of higher order coef- Fig. 3. Correlation plot for the coefficients, 
flcients, C30, C40, C50 and Cgo for the min- Ci and C2. 
imum biased HIJING events. 



4.4. Probability function: An estimate from simulation 

A more quantitative way to analyze fluctuations is to estimate the production of 
anticentauro events from simulation. For this, we define a probability function P, 
which is the fraction of events which have the magnitude of first order coefficient, Cq, 
greater than the threshold set by the sum of mean first order coefficient, obtained 
from HIJING events and five times the dispersion (< Co > +5a). For all the 
anticentauro events, the absolute value of the magnitude of first order coefficient is 
much greater than the set threshold but for the HIJING events we do not observe 
any event which satisfies the above criteria. As an example, for the azimuthal bin, 
= 27r/6, all the artificially generated anticentauro events have magnitude of the 
first order coefficient (Co = 42.7) which is much higher than the threshold value 
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(17.21). The corresponding statistical probability obtained from HIJING events is 

P(Co > 17.21) = (6) 

For the semi-central events the probability is 

P(Co > 16.61) = 2.2x10-3 (7) 

Similar behavior is observed for the other azimuthal bins. Thus, for the generated 
exotic events, high value of neutral pion fraction in the 77 -4> phase space is clearly 
exhibited by high magnitude of the Chebyshev coefficient whereas statistically the 
probability of formation of exotic event is either or very small. 

5. Summary 

In this study, we have described a method of Chebyshev polynomials to study the 
fluctuations in the neutral pion fraction of artificially generated anticentauro events 
for Pb+Pb collisions at LHC energies. This method is proved to be very efficient 
to measure the strength of fluctuations. Our analysis shows that the magnitude of 
higher range first order coefficients for exotic events is approximately 3 times greater 
than those of normal events. This trend is seen for other higher range coefficients. 
The detailed analysis is carried out by estimating the probability of formation of 
exotic events from the normal events. 
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